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We perform a numerical analysis of counterflow quantum turbulence of superfluid '^He with nonuni¬ 
form flows by using the vortex filament model. In recent visualization experiments nonuniform 
laminar flows of the normal fluid, namely, Hagen-Poiseuille flow and tail-flattened flow, have been 
observed. Tail-flattened flow is a novel laminar flow in which the outer part of the Hagen-Poiseuille 
flow becomes fiat. In our simulation, the velocity field of the normal fluid is prescribed to be two 
nonuniform profiles. This work addresses a square channel to obtain important physics not revealed 
in the preceding numerical works. In the studies of the two profiles we analyze the statistics of the 
physical quantities. Under Hagen-Poiseuille flow, inhomogeneous quantum turbulence appears as 
a statistically steady state. The vortex tangle shows a characteristic space-time oscillation. Under 
tail-flattened flow, the nature of the quantum turbulence depends strongly on that flatness. Vortex 
line density increases significantly as the profile becomes flatter, being saturated above a certain 
flatness. The inhomogeneity is significantly reduced in comparison to the case of Hagen-Poiseuille 
flow. Investigating the behavior of quantized vortices reveals that tail-flattened flow is an intermedi¬ 
ate state between Hagen-Poiseuille flow and uniform flow. In both profiles we obtain a characteristic 
inhomogeneity in the physical quantities, which suggests that a boundary layer of superfluid appears 
near a solid boundary. The vortex tangle produces a velocity field opposite to the applied superfiuid 
flow, and, consequently, the superfluid flow becomes smaller than the applied one. 

PACS numbers: 67.25.dk, 67.25.dm 


I. INTRODUCTION 

Quantum turbulence in counterflow of He II has left 
a mystery: the transition between two different kinds 
of superfluid turbulence, namely, the TI-T2 transitioni. 
Recent developments in visualization techniques have en¬ 
abled us to study this mystery. However, our theoretical 
understanding has not yet been developed enough for two 
chief reasons: 1. the difficulty in taking proper account 
of the inhomogeneity of the system and the geometry of 
the channel in the formulation and 2. the need to solve 
the coupled dynamics of the two fluids to understand the 
phenomena fully. In this paper, in order to address the 
first reason we perform a numerical analysis of the vortex 
dynamics by prescribing realistic nonuniform profiles to 
the flow of the normal fluid. The present work reveals 
characteristic inhomogeneous structure and dynamics of 
the vortex tangle that has not been reported in preceding 
works with a uniform flow profile^ We do not address 
directly the T1-T2 transition. If we can address the cou¬ 
pled dynamics of the two fluids, we will be able to reveal 
the T1-T2 transition, and we are going to report on this 
work in the near future. 

This study is motivated by the recent observation by 
Marakov et al^ of the flow profiles of the normal fluid 
in thermal counterflow in a square channel. They have 
observed the flow profiles of the normal fluid by follow¬ 
ing the motion of seeded metastable He^ molecules by 
a laser-induced-fluoresence technique^. Then in a ther¬ 
mal counterflow they prepare thin He^ molecular tracer 


lines perpendicular to the flow created by femtosecond- 
laser field ionization. Above 1 K the molecules inter¬ 
act strongly with the normal fluid and trace the mo¬ 
tion. Thus how the tracer lines move reveals the flow 
profile of the normal fluid. First, they presented evi¬ 
dence that the flow of the normal fluid is indeed turbu¬ 
lent at relatively large velocities^. In their recent stud¬ 
ies they observe the novel transition of the flow profile 
of the normal fluid^. At relatively small velocities, the 
flow of the normal fluid remains laminar, whereas the su¬ 
perfiuid is already turbulent. For heat current q < qd, 
where the value of qd obtained from the experiment is 
~50 mW/cm^, an initially straight tracer line deforms to 
a nearly parabolic shape, indicating a laminar Poiseuille 
velocity profile. For qd < q < qc 2 , where the value of qc 2 
obtained from the experiment is 80 mW/cm^, the outer 
part of the tracer line becomes flattened, which is called 
the tail-flattened velocity profile. Such a laminar profile 
has never been reported in classical fluid dynamics, and 
it must be characteristic of He H, namely, a result of the 
interaction between the normal fluid and the superfluid. 
In this work we study a numerical simulation of the vor¬ 
tex dynamics under these laminar profiles of the normal 
fluid. 

Liquid ^He shows superfluidity below a A point of 
2.17 K. In a two-fluid model one considers superfluid 
^He as an intimate mixture of normal fluid and superfluid 
components: the normal fluid has a viscosity, whereas the 
superfluid is inviscid. The velocities of the normal fluid 
and the superfluid component are denoted by Vn and Vg, 
respectively, and the densities are denoted by pn and ps, 
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respectively. The total density p = + ps is almost 

independent of temperature. The ratio of the two fluids 
depends on temperature, and Ps/p increases as temper¬ 
ature decreases. The velocity fields obey different two 
equations of motion^, and they are basically decoupled. 

The idea of quantized circulation in superfluid helium 
was suggested by Onsager— and has been observed ex¬ 
perimentally by VinenS. In a superfluid any rotational 
motion is sustained only by quantized vortices, which 
have the quantum circulation k = h/m 4 , where h is 
Planck’s constant and 77 x 4 is the mass of a "^He atom. 
The elementary excitations forming the normal fluid are 
strongly scattered by these quantized vortices. Thus, if 
there is a relative velocity between the normal fluid and 
the quantized vortices, a frictional force works between 
them, being called a mutual friction force. In the super¬ 
fluid turbulent state the mutual friction term is added to 
the equations of motion of the two fluids, so that the two 
fluids become coupled. 

Thermal counterflow in He II, which is internal con¬ 
vection of the two fluids, has been studied well as a pro¬ 
totype of quantum turbulence. One end of the channel is 
connected with a He H bath and the other end is closed. 
Upon heating the closed end, the normal fluid flows to 
the cooler side to transfer heat. The superfluid flows op¬ 
positely to the warmer end to satisfy mass conservation: 

J {PnVn + PsVs)dS = 0, (1) 

where the integral is performed over the cross section 
of the channel. Thus a relative velocity Vns = I'l’n ~ 'I’sl 
occurs between the two fluids, where the overline denotes 
the spatial average over the channel cross section. When 
the counterflow velocity exceeds a critical value, a self- 
sustaining tangle of quantized vortices appears to form 
superfluid turbulence. The measurementsiS show that 
the vortex line density L follows the relation 

= "f{Vns - Vo), ( 2 ) 

where 7 is a parameter depending on temperature and 
vq is a small parameter. 

The scheme for understanding quantum turbulence 
from the vortex dynamics was addressed by Vinen. By 
assuming homogeneous turbulence he estimated the vor¬ 
tex growth by using a dimensional analysis and modeled 
the decay process phenomenologicallyiS. He showed that 
the dynamics of the vortex tangle is described by Vinen’s 
equation 

dlj ^2 / \ 

= XlVnsL^ - X 2 L , (3) 

where xi and X 2 are temperature-dependent parameters. 
If the vortex tangle is in a steady state, dL/dt = 0, re¬ 
sulting in Eq. @. 

Schwarz has investigated counterflow quantum turbu¬ 
lence using the vortex filament model^. The observable 


quantities obtained in his calculation agree with the ex¬ 
perimental results for the steady state of vortex tangles. 
However, his simulation could sustain the steady vortex 
tangle only through some artificial procedure. Adachi et 
al. argued that this comes from the localized induction 
approximation (LIA) in which the interaction between 
vortices is neglected^. By performing the full Biot-Savart 
calculation they removed the difficulty and successfully 
obtained the steady state consistent with the experimen¬ 
tal results. We note that in these numerical works the 
system was assumed to be uniform and the physics com¬ 
ing from the inhomogeneity was missed; i.e., the effects 
of the channel walls and the resulting flow profile were 
missed. 

Counterflow quantum turbulence is known to depend 
on the aspect ratio of the cross section of the channel^. In 
low-aspect-ratio channels, the system has two turbulent 
states. An increase in the counterflow velocity is observed 
to change the laminar state to the first turbulent state 
Tl, and next to the second turbulent state T2. Melotte 
and Barenghiii suggested that the transition from Tl to 
T2 is caused by the transition of the normal fluid from 
laminar to turbulent. They addressed the issue of the sta¬ 
bility of the normal fluid with the relation to the T1-T2 
transition. In high-aspect-ratio channels, the counterflow 
exhibits only a single turbulent state T3, but there has 
been little information about what T3 is. 

The recent visualization experiments have made im¬ 
portant breakthroughs for these issues. By using micron¬ 
sized solid hydrogen tracers the Maryland groupd^ suc¬ 
ceeded in visualizing quantized vortices in superfluid “^He. 
Then they applied the technique to thermal counterflow 
to confirm directly the behavior of the two-fluid modeU^. 
Zhang et al. have studied a counterflow channel contain¬ 
ing a cylinder and showed that macroscopic eddies appear 
both downstream and upstream of the cylinder—. La 
Mantia et al. have investigated the Lagrangian dynamics 
of micrometer-sized solid particles in thermal counterflow 
and obtained the normalized probability functions of the 
particle velocities and accelerations^^. Recently, visual¬ 
ization of normal fluid flow in thermal counterflow has 
been performed. Marakov et alA found that the nor¬ 
mal fluid in a thermal counterflow has two characteristic 
laminar flows, as described in the second paragraph of 
this introduction. Recent developments in cryogenic vi¬ 
sualization techniques have shed light on the problems of 
inhomogeneous quantum turbulence. 

We mention previous numerical studies of thermal 
counterflow in He H with a nonuniform profile. Samuels^^ 
studied the dynamics of quantized vortices in laminar 
circular pipe flow with normal fluid Poiseuille flow and 
superfluid flat flow in the same direction. Aarts et alJ^ 
studied the vortex tangle in thermal counterflow by as¬ 
suming a normal-component velocity profile ether flat 
or parabolic in variously shaped channels with the LIA. 
Baggaley et al— studied thermal counterflow between 
two parallel plates by using the prescribed Poiseuille and 
turbulent profiles for the normal fluid flow. 



3 


A vortex tangle obtained by Baggaley et alr^ concen¬ 
trated near the channel walls and showed characteris¬ 
tic inhomogeneity. The simulation generated quantum 
turbulence under laminar Poiseuille flow and turbulent 
flow of the normal fluid. By examining results, they ar¬ 
gued that their results supported the scenario proposed 
by Melotte and Barenghi. A better understanding of T1 
and T2 states would be obtained by studying the flow in 
a low-aspect-ratio channel where all boundaries are solid 
except for the flow direction. This is because the T1 
and T2 states are actually observed in low-aspect-ratio 
channels and another turbulent state T3 is observed in 
high-aspect-ratio channels^. 

In this paper we performed a numerical analysis of 
thermal counterflow in a square channel. It is neces¬ 
sary to address the coupled dynamics of the two fluids 
to understand the T1-T2 transition. In a typical model 
of the coupled dynamics one would suppose that the su¬ 
perfluid is described by the vortex filament model, the 
normal fluid is described by the Navier-Stokes equation, 
and they are coupled through mutual friction^S. How¬ 
ever, it is difficult to solve fully the coupled dynam¬ 
ics. As the first essential step, we address the TI state 
where the normal fluid must be laminar. The prelim¬ 
inary results for Hagen-Poiseuille flow was reported in 
the previous study^l. The normal fluid is prescribed to 
be two types of realistic flow: Hagen-Poiseuille flow and 
tail-flattened flow. In both profiles, an inhomogeneous 
vortex tangle is obtained as a statistically steady state, 
where a boundary layer of the superfluid appears near 
the channel walls. Under Hagen-Poiseuille flow the vor¬ 
tex tangle shows a characteristic space-time oscillation. 
Tail-flattened flow was just observed by Marakov et alA, 
and its nature is unknown. The present simulation re¬ 
veals that tail-flattened flow is some intermediate state 
between Hagen-Poiseuille flow and uniform flow. 

The contents of this paper are as follows. Section H 
describes the formulation of this work. In Sec. HI we 
show results of our simulation with Hagen-Poiseuille flow. 
In Sec. IV we show results of our simulation with tail- 
flattened flow. Section V is devoted to a discussion. Sec¬ 
tion VI presents our conclusions. 


II. FORMULATION 


Here we describe the formulation of a vortex filament 
model and its numerical analysis. The equation of motion 
of a vortex filament is first explained, then the velocity 
fields of the normal fluid are described and the condi¬ 
tions of our simulation are noted. The physical quan¬ 
tities necessary for characterizing the vortex tangle are 
also introduced. 


A. Equation of motion 


In a vortex filament model^^ a quantized vortex is rep¬ 
resented by a filament passing through a fluid and has a 
definite vorticity. This approximation is very suitable in 
He H, because the core size of a quantized vortex is much 
smaller than any other characteristic length scale. 

At zero temperature the vortex filament moves with 
the superfluid velocity 


Vs = Vs,u, + Vs,b + Vs,a, 


( 4 ) 


where Vs,ui is the velocity field produced by vortex fil¬ 
aments, Vs,b is that produced by solid boundaries, and 
Vs,a is the applied uniform flow of the superfluid. The 
velocity field is given by the Biot-Savart law 


Vs,uiir) 


HI f {si — r) X dsi 


( 5 ) 


The filament is represented in parametric form as s = 
where Si refers to a point on the filament and 
the integration is performed along the filament. This 
work addresses the full Biot-Savart integral^. The veloc¬ 
ity field Us,6 is obtained by a simple procedure; it is just 
the field produced by an image vortex that is constructed 
by reflecting the filament into the surface and reversing 
its direction. According to Eq. ([1} we set Vs,a = Pns^ 
in thermal counterflow, where pns = Pn/Ps- 

By taking into account mutual friction, the velocity of 
a point s on the filament is given by 


s = Vg+ as' X (u„ - Vs) - a's' x [s' x (u„ - Us)], ( 6 ) 


where a and a' are the temperature-dependent coeffi¬ 
cients, and the prime denotes derivatives of s with re¬ 
spect to the coordinate ^ along the filament. 


B. Flow profiles of the normal fluid 


The velocity u„ of the normal fluid is prescribed to 
be two profiles, namely, the Hagen-Poiseuille profile Up 
and the tail-flattened ut profile. We suppose that the 
flow direction is along the x coordinate and that y, z are 
coordinates are normal to x, namely, u„ = u{y, z)x; the 
center of the channel is (y, z) = (0,0). 

In a rectangular channel the Hagen-Poiseuille profile is 
given by 


Up{y,z)=uo (- 1 )"*"' 

m—1,3,5,••• 

cosh(m7rz/2a) 


I - 


cosh(m7r6/2a) J m' 


cos{rmry/2a) 


( 7 ) 


where uq is a normalization factor and a and b are 
halves of the channel width along the y and z axes, 
respectively^^. 
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FIG. 1. (Color online) Tail-flattened profile Ut, normalized 
by the mean value, (a) Profile with h = 0.7. The color bar 
shows the value of Ut- (b) Profiles along 2 = 0 mm with 
various values of h. 

We realize the tail-flattened profile by composing the 
expression 

ut{y,z) =uoma.yi[up{y,z),hup{0,0)], (8) 

where max[A, B] refers to the larger value of the two 
arguments, and 0 < /i < 1 is a fitting parameter, which 
determines the flattened area. Figure[TKa) shows the pro¬ 
file of tail-flattened flow over the channel cross section, 
and Fig. [ijb) shows the profiles along z = 0 mm with 
various h values. Ut is equal to Up for /i = 0 , and ut be¬ 
comes uniform for h = 1, corresponding to simulation by 
Adachi et al^. The increase in h makes Ut more uniform, 
as shown in Fig. [TJb). 

C. Conditions 

The simulations are performed under the following 
conditions. We discretize the vortex lines into a num¬ 
ber of points held at a minimum space resolution of 
= 8.0 X 10“^ cm. Integration in time is achieved 
using a fourth-order Runge-Kutta scheme with time res¬ 
olution At = 1.0 X 10“'* s. The computing box is 
0.1x0.1x0.1 cm^. Periodic boundary conditions are used 


along the flow direction x, whereas solid boundary con¬ 
ditions are applied to the channel walls at y = ±0.05 cm 
and z = ±0.05 cm. The solid boundaries are assumed 
to be smooth completely, thus neither vortex pinning 
nor surface drag occurs. We reconnect two vortices ar¬ 
tificially, when the vortices approach each other more 
closely than A^. We eliminate vortices that are shorter 
than 5 x A^ = 2.4 x 10“^ cm. Temperatures^^ are 
T = 1.9 K (a = 0.21, a' = 0.009), 1.6 K (a = 0.098, a' = 
0.016), and 1.3 K (a = 0.036,a' = 0.014). The density 
of the normal fluid^i is Pn!P = 0.43 at 1.9 K, 0.16 at 
1.6 K, 0.05 at 1.3 K. The initial state consists of eight 
randomly oriented vortex rings of radius = 0.023 cm. 

D. Physical quantities 

In this subsection we introduce the physical quantities 
used in this paper. 

The vortex line density is defined as 

where the integral is performed along all vortices in the 
sample volume 17. 

The anisotropy of the vortex tangle is represented by 
dimensionless parameter—: 

^^[1 - («'( 10 ) 

Here fp represents the unit vector parallel to the flow 
direction. If the vortex tangle is completely isotropic, 
7 = 2/3. If the tangle consists entirely of curves lying in 
the plane normal to the flow direction, 7 = 1. 

To analyze the inhomogeneity of the vortex tangle, we 
divide the computational box into subvolumes by using 
a uniform 31^ Cartesian mesh. The vortex line density 
is defined at a subvolume as the local vortex line density 

where the integral is performed along all vortices in the 
subvolume 17'. The local anisotropic parameters are de¬ 
fined as 

These local quantities are always averaged over the flow 
direction x. Thus these are functions of position y, z 
normal to the flow direction, namely, L'{y, z) and 7'(y, z). 

III. HAGEN-POISEUILLE FLOW 

In this section we present results of our simulation un¬ 
der the Hagen-Poiseuille flow expressed by Eq. ([T]). We 
analyze physical quantities to discuss the dynamical and 
statistical properties of the quantum turbulence. 
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A. Dynamics 

We discuss the time development of L, and we 
show that the inhomogeneous vortex tangle is sus¬ 
tained by space-time oscillation. This is characteristic 
of the present case, which cannot appear in uniform 
counterflow^. 

Figure [2] shows the time evolution of L. The tempera¬ 
ture is T = 1.9 K, and the mean velocity of the normal 
fluid is Un = 0.6 and 0.7 cm/s. The vortex tangle devel¬ 
ops to the statistically steady state, even if the counter- 
flow is nonuniform. Fluctuations in L are much larger 
than those in uniform counterflow^. The amplitude of 
the fluctuation increases with Vn- The large fluctuation 
would imply that this system is unstable. To investigate 
the origin of the large fluctuation, we focus on the con¬ 
figuration of the vortex tangle. The large fluctuation of 
L is attributable to a space-time oscillation coming from 
the inhomogeneous dynamics of the vortex tangle. In the 
following the details of the oscillation are explained. 

Figure |3] shows the space-time pattern of the vortex 
tangle at T = 1.9 K and = 0.7 cm/s. The period 
of the space-time oscillation is about 0.7 s, and the os¬ 
cillation consists of four stages (a)-(d). In Fig. |3l[a), 
corresponding to a local minimum of L, the vortices are 
dilute, and the vortices remain only near the solid walls. 
Then the vortices near the walls invade the central region 
in Fig. ISl^b). These vortices make numerous reconnec¬ 
tions in the central region subject to the large counterflow 
in Fig. EIc). Hence L increases significantly to a local 
maximum. Eventually, in Fig. |3l[d), the Hagen-Poiseuille 
flow excludes the vortex tangle from the central region 
toward the solid walls. Thus the vortices in the central 
region become dilute, and they are absorbed by the solid 
boundaries. Even if the vortex tangles of (b) and (d) 
have the almost same values of L, the vortex configura¬ 
tions are quite different. Then the vortex tangle returns 
to the stage of Fig. ETa). The vortex tangle repeats 
the periodic motion, resulting in the large oscillation of 
Fig. [2] The vortex tangle sustained by this mechanism is 
more dilute than that in the case of uniform counterflow 



FIG. 2. Vortex line density as a function of t under Hagen- 
Poiseuille flow with T — 1.9 K. 


(a) t=8.80 s 


(b) t=8.99 s 




(c) t=9.20 s 


(d) t=9.34 s 




FIG. 3. Vortex tangles viewed along the flow direction in the 
statistically steady state for T = 1.9 K and = 0.7 cm/s, 
corresponding to the results of Fig. [2l These show one cy¬ 
cle of the large oscillation, while snapshots from (a) to (d) 
correspond a local minimum {t = 8.80 s), the middle of the 
increase (t = 8.99 s ), a local maximum {t = 9.20 s), and the 
middle of the decrease {t = 9.34 s), respectively. 


for the same T and Vns- 


B. Statistics 

We investigate the statistics of the quantum turbu¬ 
lence in a statistically steady state. In this subsection 
every physical quantity is averaged temporally over the 
statistically steady states and spatially over the compu¬ 
tational box. Typical experiments using second sound at¬ 
tenuation have observed values integrated over the chan¬ 
nel. The statistical values have been extensively investi¬ 
gated both experimentally and numerically. We discuss 
the properties of Hagen-Poiseuille normal fluid flow by 
comparing the statistical values with those obtained in 
previous studies of uniform counterflow. 

The statistically steady state is known to satisfy the 
relation of Eq. Here we regard the counterflow ve¬ 
locity Vns as the spatially averaged amplitude of Vns = 
Vn — Vg^a- Figure S] shows that L almost satisfies the 
relation. The vortex line density increases with T, be¬ 
cause the mutual friction becomes strong at higher T to 
develop quantized vortices. The critical velocity Vc be¬ 
low which the vortex tangle cannot be sustained is not 
investigated in this work, since the attempt to determine 
Vc correctly is time consuming and is not the main pur¬ 
pose of this work. Table U shows a comparison of 7 be¬ 
tween the present work, 7 hp, and the simulation, 7 uni, 
subject to uniform counterflow under a periodic bound- 
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Vns (cm $■■') 


FIG. 4. (Color online) Vortex line density averaged over 
the statistically steady state under Hagen-Poiseuille flow as a 
function of Vna- The dashed lines refer to the fitting of the 
data for each temperature. 


TABLE I. Line density coefficients 7 , fitting parameter vo and 
anisotropic parameter I. Given are numerical results 7 hp, Ihp 
and Vo obtained in this work under Hagen-Poiseuille flow and 
7 uni and /uni obtained by Adachi et which is a simulation 
subject to uniform counterflow under a periodic boundary 
condition. 


T 7hp 7uni Vo /hp /uni 

(K) (s/cm^) (s/cm^) (cm/s) 


1.3 

31 

53.5 

- 0.2 

0.765 

0.738 

1.6 

47 

109.6 

0.0 

0.814 

0.771 

1.9 

103 

140.1 

0.4 

0.878 

0.820 


ary condition^. The values of 7 uni quantitatively agree 
with those of the typical experiments^; thus we can re¬ 
gard 7 uni as a standard. The values of 7 hp are lower than 
those of 7 uni- This suggests that the Hagen-Poiseuille 
flow tends to suppress the growth of the vortex tangle. 
The fitting parameters vq of the present simulation are 
listed in Table HI These results are different from the value 
vq 0.1 cm/s obtained by the simulation under uniform 
periodic counterflowS. 

Figure [S] shows the anisotropic parameter / as a func¬ 
tion of Vns and T. The anisotropy is almost independent 
of Vns, in agreement with experiments^^. The dotted 
lines are the values of / averaged over Vns for each tem¬ 
perature. The anisotropy becomes large with T, because 
the increase in T intensifies the expansion of vortices to¬ 
ward the direction normal to the flow by mutual friction. 
In Table U we compare our results with the mean values 
of / obtained by Adachi et al^, which is a simulation 
subject to uniform counterflow in a periodic cube. The 
results with Hagen-Poiseuille flow are higher than those 
with uniform flow. Thus the Hagen-Poiseuille flow tends 
to increases the anisotropy. The difference between the 
results of this work and those from uniform counterflow 
increases with T. This is because the mutual friction 
becomes stronger with T to enhance the effect of the 
Hagen-Poiseuille profile. 


0.95 
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T=1.9 K 
1.6 K 
1.3 K 

• 

▲ 

T 

• • • • 

0.8 

_ ^ - 


0.75 


▼ ^ 

0.7 


- 

0.65 

_1_ 

J_1_1_ 1 _ 


0.8 0.9 1 1.1 1.2 1.3 1.4 


Vns (cm s'^) 

FIG. 5. (Color online) Anisotropic parameter / averaged 
over the statistically steady state under Hagen-Poiseuille flow 
as a function of Vns- The dotted lines are the values of I 
averaged over v„s for each temperature. 


C. Inhomogeneity 

We have discussed the physical quantities averaged 
over the whole volume. Here we analyze local physi¬ 
cal quantities introduced in Sec. Ill PI to investigate the 
structure of the vortex tangle. In this subsection every 
physical quantity is averaged temporally over the statis¬ 
tically steady states and spatially over the flow direction. 

Figure El shows the spatially dependence of the local 
vortex line density L' normalized by L. These are pro¬ 
files along z = 0 mm at three temperatures, T = 1.9, 1.6, 
and 1.3 K (from left to right). Within each subfigure, we 
plot the profiles for four different values of Vn, and the 
proHle is almost independent of Vn- One can see that the 
vortices concentrate near the channel walls. The region 
with large values of L' can be called a boundary layer of 
the superfluid. Such a boundary layer of the superfluid 
was already found by Baggaley et al. in their simulation 
of counterflow between parallel platesi^d^. This comes 
from analogy with the boundary layer of a viscous fluid, 
which is the layer near the surface where the effects of 
viscosity are significant and the intensity of vorticity is 
strong. As shown in Sec. IHID1 indeed the vorticity of 
the superfluid increases in the boundary layer. The vor¬ 
tex line density decreases very close the channel walls. 
This is a consequence of the quantized vortices enter¬ 
ing the surface normally to satisfy the solid boundary 
condition^^. This effect tends to exclude vortices parallel 
to the surface, leading to the decrease in L'. 

We discuss the dependence of the L' proHle on T in 
terms of the difference between the maximum and mini¬ 
mum values of L'. The nonuniform Vns makes the vortex 
tangle inhomogeneous through the mutual friction terms 
of a and a' in Eq. ([H]). Both the counterflow velocity 
profile and the mutual friction coefficients depend on T. 
First, the nonuniformity of Vns decreases with T, reduc¬ 
ing the inhomogeneity of the vortex tangle. According 
to Eq. ID), the counterflow velocity can be written as 

Vns = Vn- Vs,a = Vn + PnsVn, where Pns = Pn/Ps- The 
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(a) T=1.9 K 


(b) T=1.6 K 


(C)T=1.3 K 



y (mm) 


y (mm) 


y (mm) 


FIG. 6. (Color online) Profiles of the vortex line density L' normalized by L along z = 0 mm under Hagen-Poiseuille flow. 
Temperatures are T = 1.9, 1.6, and 1.3 K from left to right. 


(a) T=1.9 K 



(b)T=1.6 K 



(c)T=1.3 K 



FIG. 7. (Color online) ProHles of the anisotropic parameter /' along z = 0 mm under Hagen-Poiseuille flow. Temperatures 
are T — 1.9, 1.6, and 1.3 K from left to right. 


value of pns increases with T, making Vns more uniform, 
because the contribution of nonuniform Vn decreases. 
Second, the coefficient a increases with T, enhancing the 
inhomogeneity of the vortex tangle. The coefficient a' is 
a more complicated function of T. These effects could 
maximize the difference between the maximum and min¬ 
imum values of V at some temperature. According to 
our results, the temperature may be close to T = 1.6 K. 

Figure [7] shows the spatial dependence of the local 
anisotropic parameter 1'. These are prohles along 2 = 
0 mm at three temperatures, T = 1.9, 1.6, and 1.3 K 
(from left to right). Within each subfigure, we plot the 
profiles for four different values of Vn, and the profile is 
almost independent of Vn- We can see that /' values near 
the channel walls are lower than those in the central re¬ 
gion. This shows that the vortices in the boundary layer 
tend to be isotropic, compared with those in the central 
region. With increasing T, the boundary layer becomes 
wider, which is also shown in Fig. [S] At the channel walls 
the anisotropy becomes significantly larger. This is be¬ 
cause the quantized vortices enter the surface normally to 
satisfy the solid boundary condition, as discussed above. 


D. Superfluid velocity field 

It is important to know how the superfluid velocity 
field Vg behaves in quantum turbulence. This issue is 


closely related to the classification of two types of turbu¬ 
lence in He II at finite temperatures. One type is that in 
which both the superfluid and the normal fluid are tur¬ 
bulent and couple strongly, where the coupled turbulence 
exhibits the Kolmogorov law for its energy spectrum2ii2^. 
The other is, like thermal counterflow, that in which the 
mean superfluid and the normal fluid velocities are not 
necessarily the same. In any case, it is important to know 
how the superfluid follows the normal fluid through mu¬ 
tual friction. Samuels reported the velocity matching of 
the two fluids in laminar circular pipe flow with normal 
fluid Poiseuille flow and superfluid flat flow in the same 
directioiii^. In our system, the two fluids tend to mimic 
each other through mutual friction, but a forced relative 
motion in the thermal counterflow prevents their velocity 
matching. Then the nature of Vs is nontrivial. 

The flow direction component of the superfluid ve¬ 
locity field Vs is plotted along z = 0 mm in Fig. [5] at 
temperatures of T = 1.9, 1.6, and 1.3 K (from left to 
right). These values are normalized hy Vs a- Within each 
subhgure, we plot the profiles for four different velocities. 
The outstanding feature through all these data is the ve¬ 
locity reduction in the central region. The reduction is 
caused by the large velocity gradient near the channel 
walls, meaning that the intensity of vorticity near the 
channel walls is stronger. This comes from the charac¬ 
teristic profiles of L' and Also we can understand the 
velocity reduction in the central region by considering 








(a) T=1.9 K 


(b) T=1.6 K 


(c) T=1.3 K 





FIG. 8. (Color online) Flow direction component Vg of the velocity fields Vs of the superfluid along 2 = 0 mm under 
Hagen-Poiseuille flow. These values are normalized by ujTemperatures are T — 1.9, 1.6, and 1.3 K from left to right. 


the Biot-Savart law. Investigating in detail the vortex 
structure of Fig. [3]finds that in the boundary layer most 
of vortices have some small radius of curvature compara¬ 
ble to the layer width. These vortices propagate against 
Vn, and thus they create a velocity field parallel to Vn 
in the central region. Hence the quantized vortices in 
the boundary layer lead to the reduction in the central 
region. The dependence of on T appears mainly as 
the difference of the reduction integrated over the cross 
section. This will be discussed in Sec. IV Cl Finally, we 
explain the dependence of on the value of Vn- The 
velocity reduction of the superfluid becomes large with 
Vn- Obviously, this comes from the increase in L. As 
shown above, L increases with Vn, even if the profile of 
L' does not depend on Vn- Thus the intensity of the vor- 
ticity in the boundary layer increases with so that the 
reduction of nf becomes larger. 

The question we have to consider is whether or not 
the velocity reduction occurs in uniform counterflow. We 
conducted another simulation under uniform counterflow 
in a periodic cube for T = 1.9 K and Vn = 0.2 cm/s. The 
obtained profile of nf is almost equal to that of Vs,a- This 
shows that the reduction in the superfluid flow cannot be 
realized, unless there are a solid boundary and a resulting 
nonuniform flow of the normal fluid. 


IV. TAIL-FLATTENED FLOW 

Our model of the tail-flattened flow of Eq. m enables 
us to study how the vortex tangle is affected when the 
flow profile of the normal fluid is changed by control¬ 
ling h. The value of ft. = 0 gives the Hagen-Poiseuille 
flow discussed in Sec. m The flow profile observed by 
Marakov et alA looks to be represented approximately 
by our model with ft = 0.7. The uniform flow profile 
studied in most numerical works corresponds to ft = 1. 
Figure [5] shows how L grows for different values of ft for 
Vn = 0.5 cm/s and T = 1.9 K. The large-amplitude os¬ 
cillation still remains for low values of ft, whereas it dis¬ 
appears as the flow profile becomes more uniform with 
increasing ft. It is interesting to find that the statisti¬ 


cally steady values of L are almost saturated for ft > 0.7, 
although the flow profile is not yet uniform. Figures [JO] 
shows snapshots of the vortex tangles for four different 
values of ft. With increasing ft, the vortex tangle be¬ 
comes dense and homogeneous. Highly curved vortices 
appear at higher ft, meaning that reconnection occurs 
frequently. Visually, the vortex tangle is almost homo¬ 
geneous above ft = 0.7, and the two configurations of 
ft = 0.7 and ft = 0.9 look similar. We confine ourselves 
in this section to the case of ft = 0.7 and T = 1.9 K, be¬ 
cause in this case the flow profile mimics that observed 
by Marakov et al. and should give a typical transitional 
state between Hagen-Poiseuille flow and uniform flow. 

First, we investigate the statistics of the quantum tur¬ 
bulence in statistically steady states. Here every physi¬ 
cal quantity is averaged temporally over the statistically 
steady states and spatially over the computational box. 
One finds that they are similar to those of uniform coun¬ 
terflow. 

Figure fTTl shows L, which almost satisfies the relation 
of Eq. ([2|). We regard the counterflow velocity Vns as 
the spatially averaged amplitude of Vns = Vn — Va,a- The 
value of 7 tf obtained from this simulation is 176 s cm“^, 
which is larger than those of 7 hp = 103 s cm“^ and 7uni = 



FIG. 9. (Color online) Vortex line density under tail- 
flattened flow with various values of ft. The temperature is 
T = 1.9 K, and the normal fluid velocity is Vn = 0.5 cm/s. 
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(b) h=0.5 



(d) h=0.9 



FIG. 10. Snapshots of the vortex tangle in tail-flattened 
flow viewed along the flow direction with various values of h. 
These correspond to t = 4.00 s of Fig. [9] 


140.1 s cm“^. Because the tail-flattened profile with h = 
0.7 is almost uniform, the origin of the difference in 7 
might just be the solid boundary. This suggests that 
the solid boundary increases 7 . However, it is unknown 
how the solid boundary increases the value of 7 . The 
parameter vq is ~ 0.3 cm/s, which is between that of 
Hagen-Poiseuille flow and that of uniform flow^. This is 
consistent with the fact that tail-flattened flow is their 
intermediate state. 

Figure [12] shows the anisotropic parameter / versus 
the counterflow velocity Vns- The anisotropy is al¬ 
most independent of Vns, in agreement with experimen¬ 
tal observations^^. The dotted line shows the values of 



Vns (cm s'') 

FIG. 12. (Color online) Anisotropic parameter I averaged 
over the statistically steady state under tail-flattened flow as 
a function of Vns- The dotted line are the value of I averaged 
over the . The red solid line is the mean value of I obtained 
by Adachi et ali^ with T = 1.9 K. 


I averaged over Vns- The red solid line shows the mean 
values of / obtained by Adachi^ et al. These results show 
slightly higher values of I than those obtained from the 
simulation subject to uniform counterflow^. This is be¬ 
cause the Hagen-Poiseuille profile around the central re¬ 
gion increases the anisotropy. 

Second, we analyze the local physical quantities. Ev¬ 
ery physical quantity is averaged temporally over the 
statistically steady states and spatially over the flow di¬ 
rection. The characteristic inhomogeneous behavior ap¬ 
pears, whereas the physical quantities averaged over the 
whole volume are similar to those of uniform counterflow. 

Figure [T3| shows the spatially dependence of L'/L. We 
plot the profiles for five different values of and the 
profile is almost independent of This shows the dis¬ 
tribution of the vortices along z = 0 mm. L'/L decreases 
in the central region and very near the channel walls. 
The decrease in the central region is caused by the re¬ 
maining Hagen-Poiseuille profile, as discussed in Sec. IHII 
The value of L' increases in the flat region \y\ > 0.25 of 


120 
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Vns (cm s'"') 

FIG. 11. (Golor online) Vortex line density averaged over 
the statistically steady state as a function of Vns under tail- 
flattened flow. The value of 7 tf is 176 s cm“^ for h — 0.7 and 
T = 1.9 K. The dashed line refers to the fitting of the data. 




FIG. 13. (Golor online) Profile of the vortex line density L' 
normalized by L along 2 = 0 mm under tail-flattened flow. 
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y (mm) 

FIG. 14. (Color online) Profiles of the anisotropic parameter 
I' along z = 0 mm under tail-flattened flow. 

ut- At the channel wall, L' decreases, as discussed in 
Sec. mil The vortex-concentrated region near the chan¬ 
nel walls is a boundary layer in the tail-flattened flow. 
The boundary layer is broader and gentler than in the 
case of Hagen-Poiseuille flow. 

Figure [14] shows the spatial dependence of the local 
anisotropic parameter /' along z = 0 mm. We plot the 
profiles for five different values of Vm and the profile is 
almost independent of . The profile is almost uniform, 
but I' increases slightly in the central region, where the 
remaining Hagen-Poiseuille profile increases anisotropy. 
At the channel walls, I' increases significantly. This 
comes from the solid boundary condition for the super¬ 
fluid. 

Finally, we study the superfluid velocity field. The flow 
direction component of the superfluid velocity field Vg 
is plotted along z = 0 mm in Fig. [131 These values are 
normalized by The values are averaged temporally 
over the statistically steady states and spatially over the 
flow direction. We plot data for T = 1.9 K and five dif¬ 
ferent velocities. The superfluid flow is reduced in the 



y (mm) 

FIG. 15. (Color online) Flow direction component of 
the superfluid velocity filed Vg along z = 0 mm under tail- 
flattened flow. These values are normalized by 


center of the channel. A velocity reduction is similarly 
observed in Hagen-Poiseuille flow, although the two pro¬ 
files are different. The velocity reduction in the central 
region becomes greater with increasing This comes 
from the fact that L increases with Vn , so that the vortic- 
ity in the boundary layer increases with In contrast, 
the velocity in the outer region \y\ > 0.25 is insensitive 
to Vn- This might be a property of the velocity in the 
boundary layer. 


V. DISCUSSION 

We discuss five important topics by using the results 
of this work. 


A. Effects of the square channel geometry 

It is important to understand how the geometry of the 
square channel affects quantum turbulence. Of course, a 
real channel is not the periodic cube used in the preced¬ 
ing studies^. Some phenomena, such as the T1-T2 transi¬ 
tion, depend on the channel geometry. Because we have 
focused on the profiles of the physical quantities along 
z = 0 mm, properties attributed to the square channel 
geometry have been missed. We investigate the profiles 
of the physical quantities over the channel cross section 
to consider the geometry effect. 

Figure Uni shows profiles of the local vortex line density 
and the superfluid velocity over the channel cross sec¬ 
tion. The quantities are averaged over the statistically 
steady states. We compare the typical results of Hagen- 
Poiseuille flow and tail-flattened flow. Figures IW a) and 
fTBf b) show the profiles of the local vortex line density 
L'/L averaged over the flow direction. Under Hagen- 
Poiseuille flow the vortices concentrate strongly in the 
channel corners. This is a consequence of the quantized 
vortices tending to remain in the corners, because mutual 
friction is weaker in the region and the vortices neither 
like to expand nor shrink in the corners. In contrast, un¬ 
der tail-flattened flow the vortex line density is almost 
uniform in the flat region. In this system, vortices cease 
to remain in the corners, because mutual friction is al¬ 
most uniform over that region. Figures 116( c) and 116( d) 
show profiles of the flow direction component of the 
superfluid velocity Vg. Under Hagen-Poiseuille flow the 
value of Vg becomes a minimum in the center of the chan¬ 
nel, and increases with approaching the channel walls. In 
the corners, the superfluid velocity becomes larger than 
the applied one. Under tail-flattened flow the velocity 
reduction becomes greater in the central region than in 
the case of Hagen-Poiseuille flow. 
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vortex line density L'/L 


Hagen-Poiseuille flow Tail-flattened flow 



superfluid velocity '^s/'^s,a 


Hagen-Poiseuille flow Tail-flattened flow 



FIG. 16. (Color online) Profiles of vortex line density L' and 
the flow direction component vj of the superfiuid velocity Vs 
over the channel cross section. The temperature is 1.9 K for 
all figures. The normal fluid velocity is v„ = 0.75 cm/s for 
Hagen-Poiseuille flow and Vn = 0.5 cm/s for tail-flattened 
flow. 
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FIG. 17. (Color online) Profiles of velocity s of quantized 
vortices. The left figures show the component Sx along the 
flow direction. The right figures plot the vector field {sy, Sz)- 
The temperature is 1.9 K for all figures. The normal fluid 
velocity is Vn = 0.75 cm/s for Hagen-Poiseuille flow and v„ = 
0.5 cm/s for tail-flattened flow. The values of (sj -I- 
averaged over the whole volume are about 0.13 cm/s for (b) 
and 0.05 cm/s for (d). 


B. Velocity of vortices 

It is also important to know how the vortices move 
under the nonuniform counterflow. Thus we calculated 
the velocity s of the vortices [Eq. ([5])] over the chan¬ 
nel cross section. Figures [13 (a) and (b) show Sx and 
{sy,Sz) of Hagen-Poiseuille flow. These data was ob¬ 
tained by the average of the statistically steady states 
and the average over the flow direction. The average of 
{Sy + over the cross section is about 0.13 cm/s. 

This figure reveals the characteristic motion of vortices. 
Since Sx is negative everywhere in Fig. [T3 (a), the 
vortices turn out to move basically along v^^a- Since 
Vs,a = —{pn/Ps)vn = —0.57 cm/s, the vortices move 
slower than Vs,a in the center but faster near the channel 
walls. As shown in Fig. |I7|(b), the vortices like to move 
from the center towards the walls. The similar behavior 
appears in Fig. [T3 (c) and (d) of tail-flattened flow too. 
However, the velocity amplitude is generally smaller than 
that of Hagen-Poiseuille flow. 


C. Modification of counterflow velocity 

In Secs. m and m we regarded the counterflow ve¬ 
locity as Vns = Vn — Vs,a in the analysis of 7 . How¬ 
ever, the actual relative velocity obtained from Eq. m is 


Vns = Vn — Vs, and, in nonuniform counterflow, TJ yf Vs,a, 
as shown in those sections. Thus we should recalculate 
the counterflow velocity by taking account of the reduc¬ 
tion in the velocity of the superfluid. 

Figure [TH| shows the flow direction component of 
the superfluid velocity Vs averaged over the channel cross 
section and over the statistically steady state. Figure 
IW a) show the case under Hagen-Poiseuille flow, and Fig. 
IW b) shows that under tail-flattened flow. In Fig. IW a). 


(a) Hagen-Poiseuille flow (b) Tail-flattened flow 
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FIG. 18. (Color online) Flow direction component wf of the 
superfiuid velocity field Vs averaged over the channel cross 
section under (a) Hagen-Poiseuille flow and (b) tail-flattened 
flow. The value of uj decreases with Vn, meaning that break¬ 
down of mass conservation becomes larger. 
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TABLE II. Modified line density coefficients 7 . Given are 
numerical results with Hagen-Poiseuille flow, 7 j^p, and with 
tail-flattened flow, 7 tf, obtained in this work. 


T 

7hp 

7tf 

(K) (s/cm^) 

(s/cm^) 

1.3 

32 

— 

1.6 

50 

— 

1.9 

123 

189 


the values of Vs/vg^a for T = 1.6 K tend to be larger 
than those for T = 1.3 and 1.9 K. The values generally 
decrease from unity, and this means that mass conserva¬ 
tion [Eq. ([IJ] is broken. The breakdown becomes larger 
as Vn increases, because Vs/vs,a decreases with Vn- This 
breakdown shows that the formulation of this work is not 
sufficient for correctly describing thermal counterflow in 
a realistic channel. Thus some procedure is required to 
dynamically modify the velocities of the two fluids to sat¬ 
isfy mass conservation. 

By using these results we calculate Vns = 

Table im lists the values of modified 7 for Hagen-Poiseuille 
flow, 7 ^p, and for tail-flattened flow, 7 *f. The values of 
7 * are larger than the values of 7 , because the reduction 
in Vg becomes greater with increasing Vn- The difference 
between 7 and 7 * becomes smaller as temperature de¬ 
creases. This is because = pgjpn increases with T, 

so that the contribution from the reduction in Vg becomes 
smaller. 


D. Aspect ratio of channel cross section 

The present work addresses only the case in which the 
aspect ratio of the channel is unity, but the change of 
the ratio should seriously affect the profile of the vortex 
line density. The aspect ratio of the channel plays an 
important role in the density and the profile of the vor¬ 
tex tangle. As long as the counterflow velocity is fixed, 
the increase of the aspect ratio from unity reduces the 
counterflow velocity gradient along the long side of the 
cross section. Thus the exclusion of the vortex tangle 
shown in Figs. [31Jd) does not occur so much along the 
long side of the cross section. Hence the vortex line den¬ 
sity should increase compared to the vortex line density 
in a low-aspect-ratio channel. 

It would be meaningful to consider how the aspect 
ratio affects the normal fluid, though it is prescribed 
in this formulation. A linear stability analysis of the 
Navier-Stokes equation shows that the critical Reynolds 
number for turbulence transition of a viscous fluid in¬ 
creases significantly with decreasing aspect ratio^^. We 
expect that the normal fluid in counterflow could remain 
laminar in a low-aspect-ratio channel even if the super¬ 
fluid becomes turbulent and the vortex tangle disturbs 
the normal fluid, which may correspond to the TI state. 
To understand the Tl state, therefore, studies in a low- 


aspect-ratio channel will be indispensable. 


E. T1-T2 transition 

Considering the report^ by Marakov et al. we would 
expect the T1-T2 transition to be a transition of the nor¬ 
mal fluid from tail-flattened flow to turbulent flow. Here 
we discuss the T1-T2 transition under this assumption. 
The nature of the Tl state can be understood by using 
the results of tail-flattened flow obtained in this work. 
The Tl state will be the about homogeneous quantum 
turbulence. Under tail-flattened flow, the vortex line den¬ 
sity increases with flat parameter h, becoming saturated 
at ft, > 0.7. This suggests that the vortex line density 
does not jump at the transition from tail-flattened flow 
to turbulent flow. According to experimentsi, at the Tl- 
T2 transition the vortex line density does not jump, in 
agreement with our results. The value of 7 tf obtained in 
this work is larger than 7 of the Tl state obtained in the 
experiment. The origin of this difference is unknown, but 
addressing this issue by using coupled equations for the 
two fluids would resolve this discrepancy. 


VI. CONCLUSIONS 

In this study, we have investigated counterflow quan¬ 
tum turbulence with nonuniform flows of the normal fluid 
using the vortex filament model. The velocity field of the 
normal fluid was prescribed to consist of two nonuniform 
profiles, and the solid boundary condition is applied to 
the walls of a square channel. 

First, we conducted a simulation with Hagen-Poiseuille 
flow. The vortex line density has a fluctuation with a 
large amplitude. We investigated the configuration of the 
vortex tangle in the oscillation, and it was revealed that 
one cycle of the oscillation consists of four characteristic 
states. The relation = 'y{vng — r’o) is almost sat¬ 
isfied, but there are quantitative discrepancies between 
the results of this work and those found in uniform coun¬ 
terflow. The anisotropy found here is larger than in the 
case of uniform counterflow. Then we investigated the in¬ 
homogeneity of the vortex tangle. The quantized vortices 
concentrate near the channel walls, and the vortex tangle 
becomes more isotropic near the solid boundaries. The 
region near the solid boundary is a boundary layer of the 
superfluid, which was also found by Baggaley et aZ.— 

We investigated the superfluid velocity field. It is found 
that a velocity field Vg^gj opposite to the applied flow Vg^a 
is created by the inhomogeneous vortex tangle. Thus the 
superfluid flow becomes smaller than the applied one. 

Second, we conducted a simulation with tail-flattened 
flow. We showed that the vortex line density increases 
significantly as the flat parameter ft increases, becoming 
saturated at ft > 0.7. We confined ourselves to ft = 0.7 
to be consistent with the experiments. The value of 7 
was found to be larger than those obtained from Simula- 
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tions with uniform counterflow, and the anisotropy was 
found to be slightly larger than in the case of uniform 
counterflow. The inhomogeneity is signihcantly reduced 
in comparison to the case of Hagen-Poiseuille flow. In 
the center of the channel the vortex line density becomes 
slightly lower, and the vortex tangle becomes slightly 
anisotropic. This is because the Hagen-Poiseuille profile 
remains in the central region of the tail-flattened profile. 
The velocity field Vg^ui is different from that in the case of 
Hagen-Poiseuille flow, where the reduction of the central 
region becomes greater. 

The present simulation revealed that tail-flattened flow 
is some intermediate state between Hagen-Poiseuille flow 
and uniform flow. When we focus on the physics quanti¬ 


ties averaged over the whole volume, we find that they are 
similar to those of uniform counterflow. However, if we 
look at the spatial dependence of these values, they show 
inhomogeneous behavior. These considerations suggest 
that some in situ observation is necessary to understand 
tail-flattened flow. 
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